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The recent discovery of high-temperature superconductivity in single-layer iron selenide has gen¬ 
erated significant experimental interest for optimizing the superconducting properties of iron-based 
superconductors through the lattice modification. For simulating the similar effect by changing the 
chemical composition due to S doping, we investigate the superconducting properties of high-quality 
single crystals of FeSei-^Sa, (*=0, 0.04, 0.09, and 0.11) using magnetization, resistivity, the London 
penetration depth, and low temperature specific heat measurements. We show that the introduc¬ 
tion of S to FeSe enhances the superconducting transition temperature T c , anisotropy, upper critical 
field H c 2 , and critical current density J c . The upper critical field H C 2 (T ) and its anisotropy are 
strongly temperature dependent, indicating a multiband superconductivity in this system. Through 
the measurements and analysis of the London penetration depth A a b{T) and specific heat, we show 
clear evidence for strong coupling two-gap s-wave superconductivity. The temperature-dependence 
of A ab(T) calculated from the lower critical field and electronic specific heat can be well described 
by using a two-band model with s-wave-like gaps. We find that a d- wave and single-gap BCS theory 
under the weak-coupling approach can not describe our experiments. The change of specific heat 
induced by the magnetic field can be understood only in terms of multiband superconductivity. 

PACS numbers: 74.25.Bt, 74.25.Dw, 74.25.Jb, 65.40.Ba 


I. INTRODUCTION 


The discovery of superconductivity with transition 
temperatures of up to 55 K in iron-based pnictides has 
been at the forefront of interest over the last few years 1,2 . 
One of the most puzzling issues for these materials is the 
symmetry of the superconducting (SC) state 3 . The pair¬ 
ing symmetry in Fe-based pnictides is under debate and 
various scenarios are being considered. Among these ma¬ 
terials, iron selenide, FeSe, is of particular interest due to 
the following reasons: (i) it is considered as the simplest 
Fe-based superconductor 4 for studying the pairing mech¬ 
anism; (ii) in this system, the static magnetism is missing 
and the structural (ss 90 K) and SC (« 10 K) transition 
temperatures are well separated from each other 5 . From 
,7 Se NMR measurements, the SC transition was found to 
correlate with the enhancement of the spin fluctuations 
at low temperatures 6 ; (iii) it is characterized by the ab¬ 
sence of nesting between the hole and electron pockets 
of the Fermi surface'; (iv) the application of pressure 
around 9 GPa leads to a strong enhancement of transi¬ 


tion temperature ( T c ) up to 37 K 8 ; (vi) in this system 
the largest SC gap has been obtained by angle-resolved 
photoemission spectroscopy, which likely closes at 70 K in 
extremely tensile strained FeSe 9 . The most mysterious 
property here is not even the pressure or strain induced 
T c increase (the cuprates have already shown the ten¬ 
dency of increased T c with reduction of the dimensional¬ 
ity) , but a giant enhancement of the superconductivity at 
the Fe/SrTiC >3 interface, where SrTiOs (STO) has noth¬ 
ing in common with magnetic interaction. It seems that 
SrTiOs provides phonons that enhance superconductiv¬ 
ity in single-layered FeSe 10 . Further transport measure¬ 
ments of the single FeSe/STO has shown zero resistance 
state onset above 100 K 11 , far above the liquid nitrogen 
boiling temperature. 

Although FeSe system possesses many attractive fea¬ 
tures, the investigation of its physical properties is still in 
infancy. The material is composed of primarily PbO-type 
tetragonal FeSei-^ ( PA/nmm ) and partly of NiAs-type 
hexagonal FeSe (P63/mmc) 12 . The tetragonal struc¬ 
ture is found to transform into an orthorhombic phase 
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at low temperatures 13 . It remains unclear which of these 
phases is a superconducting one. It should be noted that 
the isotope effect experiments in Fe-based superconduc¬ 
tors 14 , show the iron isotope exponent (a c ) values be¬ 
tween 0.35 up to 0.4. Thus, one could infer that electron 
pairing in superconductors of the FeSe family is facili¬ 
tated by electron-phonon interaction. Furthermore, pure 
magnetic or spin-orbital interactions affect the interband 
coupling leading to decrease of thermodynamic T c like in 
the case of MgB 2 . Identifying the origin of the SC pair¬ 
ing mechanism is the key to understanding these inter¬ 
esting properties of FeSe. There is no general consensus 
regarding the nature of pairing at the moment. For in¬ 
stance, the bulk probes, such as specific heat 15 , Andreev 
reflections spectroscopy 16 , thermal conductivity 17 , and 
the London penetration depth A^(T 1 ) 18-20 point to the 
existence of two-gap nodeless superconductivity. On the 
other hand, evidence for nodal superconductivity in FeSe 
has been reported from the surface-sensitive probes, such 
as scanning tunneling spectroscopy 21 . 

In this paper we report on magnetization, resistiv¬ 
ity, London penetration depth, and low-temperature spe¬ 
cific heat measurements of FeSei-^Sa,. Although, simi¬ 
lar investigations have been performed in detail on anal¬ 
ogous compounds, i.e. Fe(Se,Te) 18,19,22 ’ 23 , such studies 
are lacking in the case of S-doped FeSe. Exploring the 
symmetry and structure of the order parameter, and the 
evolution of the SC gap with S doping in FeSei-^S^ sys¬ 
tem through systematic measurements of temperature 
dependent specific heat and SC penetration depth is thus 
highly desired. In order to better understand the SC 
pairing mechanism it is necessary to examine how these 
properties are affected by a different chemical composi¬ 
tion. In the first part, we deal with magnetic measure¬ 
ments in magnetic dc fields applied parallel to the c axis. 
We show that the introduction of S to FeSe enhances the 
upper critical field H c 2 , critical current density J c , and 
the T c . The upper critical fields H c 2 for H\\c and H\\ab 
have been determined from our detailed AC magnetiza¬ 
tion and specific heat studies, yielding a high supercon¬ 
ducting anisotropy T ~ 4 for x = 0.04. The anisotropy 
r of the critical field is largest close to T c and decreases 
with decreasing temperature. From the measured tem¬ 
perature dependence of the specific heat, reliable values 
of the normal-state Sonnnerfeld coefficients are obtained 
for these materials. The second part of the paper is de¬ 
voted to the study of the currently debated issue of the 
SC pairing symmetry by using high-quality single crys¬ 
tals of FeSei-^Sa,. Based on the comprehensive low-T 
measurements of the magnetic penetration depth and 
specific heat, we provide evidence for strongly-coupled 
multiband and nodeless superconductivity in FeSe fam¬ 
ily. In addition, the presence of multiple kinks in X~ b (T) 
gives strong evidence for existence of two energy gaps in 
Fe(Se,S), which implies that several sheets of the Fermi 
surface contribute to the formation of Cooper pairs. Al¬ 
though the electronic specific heat of the SC state can be 
well described by using a two-band model, the change of 


specific heat induced by a magnetic field can be under¬ 
stood only in terms of multiband superconductivity. 


II. EXPERIMENTAL DETAILS 

All preparation steps like weighing, mixing, grinding 
and storage were carried out in an Ar-filled glove-box 
(0 2 and H 2 0 level less than 0.1 ppm). FeSei-^S^ single 
crystals were grown in evacuated quartz ampoules using 
the AICI 3 /KCI flux technique in a temperature gradi¬ 
ent (a hot part of the ampule at 400°C and s cold part 
at 350°C) for 45 days 16 . The chemical composition of 
crystals was studied with a digital scanning electron mi¬ 
croscope TESCAN Vega II XMU 16 . The standard de¬ 
viation of the average S concentration allows to judge 
upon the homogeneity of S within the crystals. There¬ 
fore, the composition and in particular the S-doping 
level was obtained by an average of over several differ¬ 
ent points of each single crystal. The analysis showed 
that the approximate chemical compositions are FeSei_^, 
Fe(Seo.96±o.oiSo.o4±o.oi)i-5j Fe(Seo.gi±o.oiSo.o9±o.oi)i-5; 
and Fe(Se 0 . 89 ±o.oiSo.n±o.oi)i- 5 - The crystals have a 
plate-like shape with the c-axis oriented perpendicular 
to the crystal plane. The crystals have only a tetragonal 
/3-FeSe phase present. The lattice parameters c = (5.52 
± 0.01)1 and a = (3.77 ± 0.01)1 are found for FeSe sin¬ 
gle crystal. The quality of the grown single crystals was 
investigated by complementary techniques. 

Magnetization measurements were performed by using 
a Quantum Design SQUID. The temperature dependent 
electronic transport was measured by using a standard 
four-probe alternating current dc method within a cur¬ 
rent applied parallel to the ab plane. Electrical contacts 
parallel to the ab plane were made using thin copper 
wires attached to the sample with silver epoxy. The low- 
T specific heat was measured in the Quantum Design’s 
Physical Property Measurement System within T range 
from 2 to 14 K in magnetic fields of up to H = 9 T applied 
along c and a 6 -axis of the crystal. During the heat capac¬ 
ity measurements, the sample was cooled to the lowest 
temperature with an applied magnetic field [field cooled 
(FC)] and the specific heat data were collected using the 
adiabatic thermal relaxation technique. 


III. RESULTS AND DISCUSSIONS 
A. Magnetization 

1. DC magnetization measurements 

In Fig. 1(a), we show the magnetic susceptibility x, 
measured with zero field cooling (ZFC) and field cool¬ 
ing (FC). T c has been determined from the onset of dia¬ 
magnetic response to be around ~ 8.5, 9.58, 10.1, and 
10.7 K for x = 0, 0.04, 0.09, and 0.11 respectively. The 
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FIG. 1: (a) presents the T-dependence of the magnetic susceptibility x i n an external held of 10 Oe applied along the c-axis. 

X has been deduced from the dc magnetization measured by following ZFC and FC protocols, (b) and (c) show the isothermal 
magnetization M vs. H loops measured at different temperatures ranging from 2 to 9 K up to 9 T applied along the c-axis. 
(d) illustrates the T-dependence of the critical current density J c values at H = 0 for x = 0, 0.04, 0.09 and 0.11. (e) shows 
the temperature dependence of the in-plane electrical resistivity in zero held up and represents a zoom of the resistivity data 
around the superconducting transition. 


FC and ZFC data show a sharp diamagnetic signal on¬ 
set. The SC volume fraction of the crystals is close to 
1, thus confirming bulk superconductivity and the high 
quality of the investigated systems. The clear irreversibil¬ 
ity between FC and ZFC measurements is consequence 
of a strong vortex trapping mechanism, either by surface 
barriers or bulk pinning. The fact that the hysteresis 
loops are symmetric around M = 0 line, points to rel¬ 
atively weak surface barriers and is indicative of strong 
bulk pinning 24 . This consideration holds for all studied 
temperatures and investigated samples, even close to T c 
and indicates that vortex penetration occurs at a field 
close to the thermodynamic H c i (corrected by the de¬ 
magnetization factor). Magnetization curves [Figs. 1(b) 
and (c)] show presence of a second peak for FeSeo. 91 So .09 
and FeSe 0 . 96 S 0.04 for H || c. The second peak effect 
has been studied extensively and its origin may be at¬ 
tributed to various mechanisms. The superconducting 
hysteresis loops M(H) exhibits no paramagnetic back¬ 
ground, which indicates that our investigated samples 
contain negligible amount of magnetic impurities and all 
Fe atoms are in nonmagnetic state of Fe 2+ . From the 
irreversibility of the magnetization hysteresis loops in 
M(H), we have extracted the magnetic field dependence 
of the critical current density J c at different temperatures 
(see Fig. 1(d)). We used Bean’s critical state model 25 in 
which the critical current is constant across the sample 
and the critical current density in a platelet sample is 
given by: 


Jc = 


20AM 

Ki-i)]’ 


(1) 


where AM = Mrf n —M„ p , Md n and M up are the magneti¬ 
zation values measured on the decreasing and increasing 
branches of M(H), respectively, a [cm] and b [cm] are 
sample sizes perpendicular to the applied field (a < b). 
Here AM is in electromagnetic units per cubic centimeter 
and the calculated J c is in Ampere per square centime¬ 
ter. The calculated J c values are summarized in Table I. 
Figure 1(e) shows the in-plane resistivity data for x = 0, 
0.04, 0.09, and 0.11 samples near T c . A sharp SC tran¬ 
sition is seen in all of the samples which is in agreement 
with the magnetization data in Fig. 1(a). Upon cool¬ 
ing the resistivity decreases monotonically and shows a 
metallic behavior. 


2. AC magnetization measurements 

Figures 2(a) and (b) depict the temperature depen¬ 
dence of the volume AC susceptibilities x'v for x = 0.04, 
and 0.09 respectively. The measurements were done in 
an AC field with an amplitude Hac = 5 Oe and a fre¬ 
quency f = 1 kHz in DC fields up to 9 T parallel to the c 
axis. Special care has been taken to correct the magne¬ 
tization data for demagnetization factor, where the de¬ 
magnetization factor has been estimated based on crys- 
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FIG. 2: (a) and (b) summarize the temperature dependence of the complex AC-susceptibility components 4ttx'v of x = 0.09, 

and 0.11 measured in an AC field with an amplitude of 5 Oe and a frequency of 1 kHz. The data were collected upon warming 
in different DC magnetic fields after cooling in a zero magnetic field, (c) presents the imaginary part of AC at various v m for 
x = 0. (d) shows the values from the position of the maxima of the imaginary part in (c) vs. i>m for x = 0, 0.04 and 0.09. 


tal dimensions 26 . In general, AC-susceptibility measure¬ 
ments can be used for an investigation of the flux dy¬ 
namics in superconductors. The imaginary part x"v is 
related to the energy dissipation in the sample due to vor¬ 
tex motion and the real part x’v is related to the amount 
of Meissner currents screening. Both functions depend 
on the ratio between the skin depth S s and the sample 
dimension L in the direction of the flux penetration. In 
the normal state S s ~ (pn//) 0 ' 5 , where p n is the normal- 
state resistivity and / is the frequency 27 . In the super¬ 
conducting state, the skin depth S s oc \l if an external 
magnetic held is below the first critical held H c \, where 
Al is the London penetration depth. For magnetic helds 
above H c i, S s oc Lp, where Lb ~ B ac /J c is the Bean’s 
penetration depth and J c is the critical-current density. 
In general, if L <C S s the AC held penetrates completely 
the sample, although the susceptibility is small. In the 
opposite case, L 5 s , most of the sample volume is 

screened. Therefore, 4n\'v = — 1 and x"v 0. 

In accordance with this, the AC-susceptibility data 
measured at low temperatures confirm the bulk super¬ 
conductivity of the investigated crystals. The transition 
temperature T c has been extracted from the bifurcation 
point between x'v and x"v This point is related to 
the change in the resistivity due to the superconduct¬ 
ing transition. It can be also used for the determina¬ 
tion of the temperature dependence of the upper critical 
held H c2 from the AC-susceptibility data measured at 
various DC helds. Therefore, the most natural way to 
investigate the vortex dynamics is to repeat y ac vs. T 
scans at a hxed Hd c at different frequencies v m in order 
to employ empirical peak functions around the maxima. 
Figure 2(c) presents the temperature dependence of the 


imaginary part of AC susceptibility at various frequen¬ 
cies v m for FeSe. One can clearly see that the peak max¬ 
ima shifts to a higher temperatures upon increasing the 
frequency which is apparently due to the motion of vor¬ 
tices. Figure 1(d) shows Tp values, the position of the 
maxima of imaginary part in (c) vs. v m . One can notice 
that, similarly to what was observed in YBa 2 Cu 30 y 28 
and CeFeAsO 0 . 92 F 0 . 08 29 ; the quantity 1 /T p presents a 
frequency dependence over the explored range of v the 
explored frequency range at hxed applied held H [H ac 
= lOe and Hd c = 1 T, see Fig. 2(d)]. The latter phe¬ 
nomenology is well verified for all the samples x = 0, 
0.04 and 0.09. The dashed line outlines the expected 
logarithmic trend typical for thermally-activated process 
according to: 1/T p (u m ) = — ^ln^. The parameter 
u 0 represents an intrabancl condensate vortex hopping. 
From the latter equation, it can be observed that the 
logarithmic behavior of 1 /T p is mainly controlled by the 
parameter Uq, playing the role of an effective depinning 
energy barrier in a thermally activated hux creep model. 

B. Specific heat 

Low-temperature specihc heat Cp , being equal to the 
temperature derivative of the entropy S, and probes the 
gap structure of bulk superconductors. The thermody¬ 
namic Cp measurements well complement the magnetic 
(A) measurements, since the former is hardly affected by 
vortex pinning. The information about the pairing sym¬ 
metry is contained in the C e i, which is proportional to 
the quasiparticle density of states (DOS) at the Fermi en¬ 
ergy. First we address the zero-held T-dependent specihc 
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FIG. 3: (a) C p /T vs. T of x = 0, 0.04, 0.09, and 0.11 in zero magnetic field. Data are shifted by an offset along the y axis 

for clarity. The upper inset: enlarged C p /T vs. T plot near the SC transition for x = 0.04. The lines illustrate how C p /T c 
and T c are estimated. The dashed line is the fitting of the specific heat below 141< by using C p = 7 n T + aT 3 + /3T 5 . (b) The 
doping dependence of the normalized specific heat jump is AC p /xnT c for x = 0, 0.04, 0.09, and 0.11. (c) and (d) summarize 
the temperature dependence of the specific heat of x = 0.04, and 0.09 respectively in various applied magnetic fields up to 
9 T parallel to the c axis and parallel to the ab plane as presented in the insets, (e) presents the temperature dependence of 
the specific heat of both orientation for 0T and 9T two data sets with the same T c value for the two directions. The inset 
highlights the electronic specific heat data after subtracting the phonon contribution for x = 0.04 of two data sets with the 
same T c value for the two directions. The data present the 0.5 T B || c and for 2T B || ab. (f) the temperature dependence of 
the complex AC-susceptibility components 4ttx'v and 4nx"v of x = 0.04 of both orientation for 0 T and 9 T measured in an AC 
field with an amplitude of 5 Oe and a frequency of 1 kHz. The data in (e) and (f) show that the electronic specific heat divided 
by temperature for B || ab and B || c and the AC-susceptibility represent an anisotropy of T = 4 for S-doping concentration. 


heat data of FeSei_a;S a . plotted as C p /T vs T (Fig. 3(a)). 
A clear sharp anomaly is observed due to the SC phase 
transition. In order to determine the specific heat re¬ 
lated to the SC phase transition we need to estimate the 
phonon (Cph) and electron (C e i) contributions. At low 
temperatures, C e i behaves linearly with temperatures, 
while C p h varies as Cph oc T 3 . In order to improve the re¬ 
liability at higher temperatures, we use a second term of 
the harmonic-lattice approximation below 14 K. The data 
can be well described by Cei + Cph = 7 T+aT 3 +/3T 5 (see 
the dashed line in Fig. 3(a)), in which the lattice contri¬ 
bution is accounted for by aT 3 + /3T 5 . The Sommerfield 
coefficient values are 5.3(1), 5.1(0.5), 4.9(0.5), and 
4.95 for x = 0, 0.04, 0.09, and 0.11 respectively. The 
estimated universal parameter A C e i/"f n T c of the specific 
heat at T c is ss 2.14, 2.43, 2.2, and 1.95 mJ/mol K 2 for 
x = 0, 0.04, 0.09, and 0.11, respectively. 

According to the BCS theory, the specific heat jump 
of a superconductor at T c should follow A C P l^ n T c = 
1.43 in the weak coupling limit. It is so far well re¬ 
ported that a reduced jump in the specific heat A C p /T c 
compared to that of a single-band s-wave superconduc¬ 
tor might be related to a pronounced multiband char¬ 
acter of superconductivity with rather different partial 
densities of states and gap values 2 ' . However, jumps 


of specific heat at T c in these materials show a devia¬ 
tion from the trend established by Bud'ko-Ni-Canfield 
(BNC) scaling A C p /T c oc T 230 . This power law seems 
to be appropriate for many iron-based superconductors, 
so far for the 122 systems, i.e., Ba(Fei_ x Co ;c ) 2 As 2 and 
Ba(Fei_ x Nij,) 2 As 2, 31 One of the possible reasons for this 
universal relation might be a strong pair breaking and the 
impurity scattering effect in a multiband superconduc- 
tor 30 ’ 32 . Recently, specific heat jump shows also a devi¬ 
ation from that trend in FeSeo.sTeo.s 33 . In addition, the 
heavily hole-doped superconductors (K,Na)Fe 2 As 2 and 
stands out from the other Fe pnictides 34 . 

From the extracted y n values we have estimated 
the values of universal parameter C e i/"f n T c in Table I, 
(Fig. 3(b)). However, a domelike dependence on dop¬ 
ing is seen similar to the one in NaFei-^Co^As 35 and 
Ba(Fei_ x Co x ) 2 As 2 36 . The values recorded in in Fig. 3(b) 
are larger than the BCS weak-coupling prediction of 
1.43. Therefore, we notice that the reduced specific 
heat jumps are larger than those of optimally doped 
Ba(Fei_ x Co :r ) 2 As 2 37 . Nevertheless, the values are com¬ 
parable to those of optimally doped Bai_ a; K x .Fe 2 As 2 38 . 
In addition, band-structure calculations in FeSe esti¬ 
mated the value of 7 0 of about 2.2 mJ/mol K 2 . Since 
7n = (1 + A ei-bosho , where A e i-bos is the total coupling 
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FIG. 4: Upper panel summarizes phase diagram of H c 2 
vs. temperature of x = 0.04 and FeSeo. 45 Teo. 55 39 (T c = 14.5 K) 
for the field applied parallel and perpendicular to c. T c has 
been estimated from an entropy-conserving construction and 
AC measurements. H c 2 in 39 is determined from transport 
measurements. The inset illustrates phase diagram of H c 2 
vs. temperature of the FeSe. The red solid points in x = 
0.04 data are estimated from the ac magnetization while the 
black data represents the specific heat for H || c and H || ab. 
Lower panel represents the phase diagram with the fit to the 
Ginzburg-Landau equation (dotted line) and the WHH model 
for A = 0, a = 0 (solid lines). The inset shows the anisotropy 
T = H^ ±c / H^ c determined by an interpolation of the H c 2 
curves. The line is a guide to the eye. 

strength between the quasiparticle (QP) and bosons 40 , 
the A e i-bos value is estimated to be « 1.4. Without any 
model fitting, the values of normalized specific heat jump 
and the A e i-bos constant have already further confirm a 
stronger electron-boson-coupling strength in FeSei-^S^. 


C. The upper critical fields H C 2(T) and their 
anisotropy 

Figures 3(c) and (d) summarize the temperature de¬ 
pendence of the specific heat measured at different ap¬ 
plied magnetic fields parallel to the c-axis (and parallel 
to the a&-plane) as shown in the insets x = 0.04 and 0.09, 
respectively. In zero-field specific-heat measurements, a 
very sharp anomaly is clearly seen. This anomaly is at¬ 


tributed to the superconducting transition at T c . This 
specific-heat jump is systematically shifted to lower tem¬ 
peratures upon applying DC magnetic fields of up to 9 T 
in both orientations. In order to determine the super¬ 
conducting transition temperature for each field, an en¬ 
tropy conserving construction has been used. The in¬ 
set in Fig. 3(e) presents the temperature dependence of 
the electronic specific heat near the transition tempera¬ 
ture at 0.5T||ab and 2T||ab data. The results yield an 
anisotropy of about T = 4 for x = 0.04. The extracted 
data at each field were used to map out the superconduct¬ 
ing phase diagram depicting H c2 (see Fig. 4). In order to 
determine the upper critical field H c2 for the c orienta¬ 
tion, we use the Ginzburg Landau (GL) equation as an 
initial step : 41 

H c2 = H c2 (0)[^}, ( 2 ) 

where t is the reduced temperature t = T/T c . The fit is 
shown with a dotted line in the inset of Fig. 3(f). Another 
possibility to extract the upper critical field H c 2 (0) is to 
consider the single-band-Werthamer-Helfand-Hohenberg 
(WHH) formula 42 for an isotropic one-band BCS super¬ 
conductor in a dirty limit. An example of WHH fit is 
shown with the solid lines in the inset of Fig. 3(f) for both 
orientations. The WHH theory (a = 0, X so = 0) predicts 
the behavior of H c2 (T c ) taking into account paramag¬ 
netic and orbital pair-breaking 42 . Here, a is the Maki 
parameter which describes the relative strength of orbital 
breaking and the limit of paramagnetism. A so (when A so 
> 0 , the effect of the spin-paramagnetic term) is the spin- 
orbit coupling strength. The orbital limiting field I?°| b 
at zero temperature is determined by by a slope at T c 
as n 0 H°^ h = 0.69 T c (d» 0 H c2 /dT) | Tc , where <J ^| Tc is 
indicated by the slopes of the fitted straight lines. It 
should be borne in mind that Fig. 4 (upper panel) re¬ 
ports in the same plot the data for FeSeo. 45 Teo. 55 39 and 
the anisotropy of H c2 (T ) is found to be 2. 

The upper critical field values at T = 0 for the FeSe 
and FeSe 0 . 96 S 0.04 have been evaluated to be ^jlqH^{0) = 
12.8(1), 17.5(1) T and /j 0 H% b) (0) = 24.4(2), 67.5(2) T. 
The anisotropy is found to be T = H^/H^ ~ 2 and 4 
for x = 0 and 0.04, respectively. The extracted values 
for H c 2 ( 0 ) are summarized in Table I for the other sam¬ 
ples. The observed small differences between the data 
obtained from the specific heat and the AC magnetiza¬ 
tion for H || c (see upper panel Fig. 4) is not surprising 
since these methods naturally imply different criteria for 
the T c determination. It is evident that, the one-band 
WHH model fails to satisfy the extracted H c2 ( 0), i.e., 
the fit shown by the red solid line in the lower panel 
of Fig. 4. Therefore, we believe that the observed de¬ 
viation from the single band WHH model is related to 
multiband effects. Furthermore, the temperature depen¬ 
dence in H c2 (T) displays a non-linear behavior and shows 
a slightly concave curvature at low temperatures. This 
behavior is reminiscent of the one reported in Fe-based 
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superconductors 2 ' ,43 in which similar H c 2 (T) curves were 
well described by a multiband effect. This claim is 
supported by the indications of a two-band-like behav¬ 
ior from the zero-field specific heat measurements and 
the London penetration depth (discussed below). The 
calculated upper critical field H c 2 (0) and average slope 
—dH^/dT values of the superconducting single-crystal 
samples are compared in Table I. 

In all of the investigated S-doped samples, both specific 
heat and AC susceptibility measurements at 9 T|| ab data 
show a sharper transition compared to the 9T data || c, 
indicating a highly anisotropic crystal, [see Figs. 3(e) and 
(f)]. From the behavior of H c 2 vs. T for the different field 
orientations we calculate the anisotropy T = c 

using a linear interpolation. The results are outlined in 
the inset of the lower panel of Fig. 4. Our data allow 
tracking of F up to temperatures very close to T c . T 
increases upon approaching T c and reaches a value of 
about 4 at the critical temperature for x = 0.04. This 
finding is in contrast to the results found in other Fe- 
based superconductors, which suggests a saturation or 
even a decrease of T close to T c 44 . This indicates that an 
orbital pair breaking is dominating the suppression of su¬ 
perconductivity close to T c . At lower temperatures, the 
isotropic Zeeman pair breaking becomes more important, 
which leads to lowering of T. Moreover, the strong In¬ 
dependent superconducting anisotropy T = H^ 2 ±c /'H^} c 
provides further evidence for multiband scenario as in 
the case of the two-band superconductor MgB 2 45 . Sur¬ 
prisingly, this anisotropy is considerably larger than the 
typical value of T ~ 2 and 2.6 found in nearly opti¬ 
mally hole-doped (BaK)Fe 2 As 2 46 and in FeSeo.45Teo.55 39 
(T c = 14.5K), but lower than the ones determined in 
SmFeAsO0.85F0.15 and La(0,F)FeAs thin films 47,48 . On 
the other hand, these values are comparable with T val¬ 
ues of e.g., KFe 2 As 2 27 and LaFePO 49 . 
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FIG. 5: The field dependence of the initial diamagnetic part of 
the magnetization curves measured at various temperatures 
for H || c in FeSei-xSa; single crystal, where x = 0, 0.04, 
0.09, respectively. The inset of (a) presents the magnetic field 
dependence of magnetization in FeSe0.91S0.09 single crystal at 
different temperatures close to T c . The insets in (b) and (c) 
depict an example used to determine the H c 1 value using the 
regression factor R, at T = 2 K. 


D. Superconducting energy-gap structure 

1. London penetration depth 

The London penetration depth A, is a fundamental 
parameter characterizing the SC condensate and probes 
the gap structure of bulk superconductors. The In¬ 
dependence of A is directly determined by the gap func¬ 
tion ALT). A(T) = A(T = 0)+5X(T) behaves as <5A(T) oc 
exp(^^) at low T reflecting a nodelsss superconducting 
gap A with s-wave symmetry. In d -wave pairing sce¬ 
nario containing line nodes, SX(T) oc T at low T in the 
clean limit. The experimental determination of the Lon¬ 
don penetration depth A(0) via measurement of the lower 
critical field H c 1 is challenging since H c 1 is an equilib¬ 
rium thermodynamic field. The temperature dependence 
of the SC penetration depth provides another method for 
detecting the existence of multiple gaps 50 . A popular ap¬ 
proach to measuring H c \ is by tracking the virgin M(H) 


curve at low fields and identifying the deviation from the 
linear Meissner response which would correspond to the 
first vortex penetration (see Fig. 5). This technique im¬ 
plicitly relies on the assumption that no surface barriers 
are present. We have confirmed the absence of the sur¬ 
face barriers in our case from the very symmetric DC 
magnetization hysteresis curves M(H) (see Fig. 1(b) and 
(c) and the inset in Fig. 5(a)). On the other hand, if sur¬ 
face barriers were predominant, the first vortex entrance 
would occur at much higher field (~ H c ). Thus absence 
of surface barriers is a very important for obtaining reli¬ 
able estimates of the thermodynamic lower critical field. 
The transition from linear to non-linear M(H ), was de¬ 
termined by a user-independent procedure consisting of 
calculating the regression coefficient I? of a linear fit to 
the data points collected between 0 and H, as a function 
of H (see the insets in Fig. 5). In contrast to tracking the 
virgin M(H) curve at low fields at several temperatures, 
in which case a heavy data post-processing is needed, (see 
Fig. 5), here a careful measurement protocol needs to be 































FIG. 6: (a-c) The normalized SC electronic specific heat of the three samples after subtracting the phonon contribution as a 

function of reduced temperature , T/T c . The inset in (b) presents the entropy in the normal and SC state as a function of T. 
(d-f): The temperature dependence of the London penetration depth for x = 0, 0.04, and 0.09, respectively. The inset of (d) 
presents the temperature dependence of the magnetic penetration depths X a b vs. T. The inset of (e) depicts an example used 
to determine the H c i value using the trapped moment at T = 2 K of the typical plot of y/Mt vs. H. The solid red line is a 
linear fit to the high-field data of y/Mt vs. H. H c \ values are determined by extrapolating the linear fit to y/Mt =0. The inset 
in (f) represents the phase diagram of H c i for the field applied along the c axis. The dashed lines represent the theoretical 
curves based on single-band weak-coupling BCS theory, while the dotted lines present the d- wave approximation. The solid 
lines represent the curves of the two s-wave gap model. 


followed with little data analysis. Indeed, the H c i values 
from the virgin magnetization data are close to those ob¬ 
tained from the onset of the trapped flux moment M t [see 
the inset in Fig. 6(e)]. Here, M t is obtained by following 
sequence: (i) warming the sample up to temperatures 
above T c i.e., 12 K, then (ii) cooling the sample in zero 
field down to the particular T, and, subsequently (iii) 
increasing the applied magnetic field to a ceratin max¬ 
imum value H m and in a last step (iv) measuring the 
remanent magnetization M t after the applied field has 
been removed. The field H m at which M t deviates from 
zero determines the H c \ value at the desired temperature. 
Then, the extrapolation y/M t —► 0 determines the exact 
value of the H c i. The inset of Fig. 6(e) presents the typi¬ 
cal plot of y/Mt vs. the applied field H, for FeSeo. 96 So .04 
single crystals. The solid line is a linear fit to the high- 
held data of y/Mt vs. H. 7J cl is determined by intercept 
of the fit with the abscissa. 

The above measured values of H c i need to be cor¬ 
rected due to the finite demagnetization effects. Indeed, 


the deflection of field lines around the sample leads to 
a more pronounced Meissner slope given by M/H a = 
— 1/(1 — N), where N is the demagnetization factor. 
Taking into account these effects, the absolute value of 
H c i can be estimated by using the relation proposed by 
Brandt 51 . For our sample we find N ss 0.96, 0.95, 0.97, 
and 0.96 for FeSeo.96So.o4, FeSe0.91S0.09, FeSeo. 89 So .11 
and FeSe respectively. In order to shed light on the pair¬ 
ing symmetry in our system, we estimated the penetra¬ 
tion depth using the traditional Ginzburg-Landau (GL) 
theory, where H cl is given by: fi 0 H^ = (</> 0 /47rA^ b ) lnK c , 
where </>o is the magnetic-flux quantum <j>o = h/e* = 
2.07 x 10 _7 Oe cm 2 , k c =X a b/^ab is the Ginzburg-Landau 
parameter. The value of k was determined from the 
equation: = - ') g 2 °‘ 5 * It should be mentioned that 

the SC penetration depth is a very important physical 
quantity and it is sensitive to the absolute value the or¬ 
der parameter(s); and, in that sense, also sensitive to any 
nodes or a deep local minima of the gap. It is worth men¬ 
tioning that the SC penetration depth is also dependent 










































9 


TABLE I: Compilation of the superconducting parameters of samples with various T c . We show the T c (K), J c (10 4 A/cm 2 ), 
A ai ,(0) ((15) nm), T = Hc 2 ±c /H ^ c , 7 „ = (mj/mol K 2 ), the universal parameter (AC e ;/ 7 nTc), — \t c (T/K), upper critical 

field H c 2 (T), /3 = 10~ 4 inJ/mol K 4 , a = 10 _ 7 mJ/mol K 6 , d-wave (Ao/fcsT c ), the superconducting gap ratio ( 71 , 72 / 771 ), and 
two s-wave gaps (an = Ai /ksT c and 02 = A 2 /ksTo) extracted for the investigated samples. 


X 

Tc 

Jo 

A of ,(0) 

r 

7« 

A Col hnT c 

dp-aHc.2 1 . 
dT \ Tc 

Ho 2 

P 

a 

d- wave 

71 , 72 / 7 ™ 

cti/ot2 (C p , \ a b) 

0 

8.5 

1.1 

446 

2 

5.3 

2.14 

2.1 

12.8 

4.34 

-0.384 

2.8 

0.4, 0.6 

0.88, 0.79 / 2.22, 2.05 

0.04 

9.58 

1.3 

372 

4 

5.1 

2.43 

2.6 

17.5 

4.8 

-3.62 

2.36 

0.44, 0.56 

1.9, 1.85 / 2.5, 2.3 

0.09 

10.1 

1.35 

433 

3.5 

4.9 

2.2 

2.7 

19 

3.6 

-2.5 

3.05 

0.47, 0.53 

2.2, 2.1 / 2.35, 2.28 

0.11 

10.7 

1.45 

415 

3.5 

4.95 

1.95 

2.73 

20.2 

4.1 

-2.9 

3.12 

0.42, 0.58 

1.96 / 2.39 


on the distribution of Fermi velocities. In this context 
it is most sensitive to fast electrons in sharp contrast to 
the upper critical field which is highly sensitive to the 
subgroup of electrons with low Fermi velocities. 

In Fig. 6(d-f), we analyze the temperature depen¬ 
dence of the London penetration depth for the sam¬ 
ples with x = 0, 0.04, and 0.09, respectively. We 
compare our data to the d -wave and single-gap BCS 
theory under the weak-coupling approach (see dotted 
and dashed lines in Fig. 6(d-f)). Indeed, both quanti¬ 
ties lead to a rather different trend and show a sys¬ 
tematic deviation from the data in the whole T-range 
below T c . On the other hand, we also apply a phe¬ 
nomenological two-gap model which is in line with the 
multigap-superconductivity reported by Carrington and 
Manzano 52 . Within this model, the temperature depen¬ 
dence of each energy gap can be approximated as: 52-54 
A i(T) = A.j(0)tanh[l.82(1.018(^/7 — l)) 051 ]. According 
to Ref. 55, for each band, A“ 2 (T) is given by: 


K\T) 


Ai(T)tanh(f^) 

A?(0)Ai(0) 


(3) 


where Aj(0) is the residual penetration depth for each 
band, k B is the Boltzmann constant. Considering dif¬ 
ferent partial contributions of each band to the over¬ 
all A(T), we use the following expression: A _2 (T) = 
rAj~ 2 (T) + (l —?’)A^ 2 (T) with r being the weighting factor 
indicating the contribution of the small gap. The best de¬ 
scription of the experimental data is obtained using val¬ 
ues of Ai /k B T c = 1.72±0.3, 1.79±0.25, and 0.79±0.15 
and A 2 /k B T c = 2.28±0.3, 2.1±0.25, and 1.95±0.2 for 
x = 0.09, 0.04, and 0, respectively. The weighting fac¬ 
tor is found to be around r = 0.25±0.08, 0.38±0.1, and 
0.22±0.2 for x = 0, 0.04, and 0.09, respectively. The 
fits are represented by solid red lines in Fig. 6(d-f). The 
extracted gap values for FeSe are comparable to those 
obtained from the two-band s-wave fit of the specific 
heat data and the Andreev reflection spectroscopy re¬ 
sults 15,16 . It is worth pointing out that the A“ & 2 (T) 
of the SC samples [see the inset of Fig. 6(d)] does not 
saturate at low temperatures, as it could be expected 


for a fully gapped clean s-wave superconductor. A a b{T) 
is nearly constant at low temperatures, which demon¬ 
strates negligible quasiparticle excitations. The above 
penetration depth results are consistent with the pres¬ 
ence of two s-wave-like gaps. Both gap values in the 
S-doped samples are considerably larger than the BCS 
weak-coupling limit. These observations show clearly 
that there are no nodes in the SC energy gap indicat¬ 
ing a strong-coupling multiband (and nodeless) super¬ 
conductivity in iron chalcogenide Fe(Se,S) superconduc¬ 
tors. A similar possible strong coupling multiband su¬ 
perconductivity in Fe(Se,Te) has been conjectured from 
a detailed penetration depth and specific heat experi¬ 
ments 22,23 . The temperature dependence of the magnetic 
penetration depth of d- wave superconducting gap calcu¬ 
lations was performed by using the following functional 
form : 18,20 

Kb( T ) _ , + I f 2n T EdEdy 

Kb (0) 71 Jo J A(T,ip) dE yjE 2 — A(T, <p) 2 ’ 

(4) 

where / = [1 + exp (E/K B T)\ 1 is the Fermi function, 
if is the angle a long the fermi surface, and A(T, <p) = 
A 0 5(T/T c )g(<p) (Ao is the maximum gap value at T=0). 
The function g(tp) is given by g d (<f) = \cos(2ip)\ for the 
d- wave gap. The results of the analysis are presented in 
Fig. 6 (d-f) by dotted lines. The fit to of the experimental 
data for the d- wave case we get for A 0 /k B T c = 2.8, 2.36, 
and 3.05 for x = 0, 0.04, and 0.09. It is obvious that 
the d- wave case cannot describe the penetration depth 
data. On the other hand, the experimental data are well 
described by the two-gap s wave models. 

It is noteworthy that in a FeTeo. 5 sSeo .42 system, a 
careful analysis of the superconducting and normal state 
properties indicates a possibility of strong coupling su¬ 
perconductivity 56 . This study is followed by precise mea¬ 
surements of the temperature dependence of the London 
penetration depth by Cho et al. Their analysis strongly 
suggest a presence of two s-wave-like gaps with magni¬ 
tudes Ai /k B T c = 1.93 and A 2 /k B T c = 0.9 57 . These 
two precise measurements 56,5 ' were followed by the com¬ 
ment of Klein et al. 5S and response of K. Cho et al. 59 . 
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In the latter case, the authors have shown convincingly 
that previous studies 22 most likely have issues with pair¬ 
breaking scattering. The authors have reported that the 
presence of strong scattering hinders any determination 
of gap values from the temperature dependence of the 
superfluid density. 


2. Specific heat 

The normalized zero-field data C e i/”/ n T as a function 
of the reduced temperature T/T Cl obtained after sub¬ 
tracting the C p h 7 is presented in Fig. 6 (a-c) together with 
the fits to various models. It is obvious from Fig. 6 (a-c) 
that the superconducting transition at T c is well pro¬ 
nounced, with a sharp jump in C e \ at T c . The entropy 
conservation required for a second-order phase transi¬ 
tion is fulfilled as shown in Fig. 6 (b). This check war¬ 
rants the thermodynamic consistency for both: the mea¬ 
sured data and the determination of C e \. We have at¬ 
tempted best fits to the data using three different mod¬ 
els: single-band weak-coupling BCS theory with the s- 
wave gap A( 0 )/IbT c = 1.76; a d -wave calculation using 
A = A( 0 ) cos(2<^>); and two-gaps s-wave in Figs. 6 (a-c). 
Below T c we observe systematic deviation of both single¬ 
gap and the d-wave fit from the data showing a higher 
jump at T c than the s-wave model. Thus we focus our 
discussion on the possibility of two SC energy gaps us¬ 
ing the generalized a-model, that explains the specific 
heat behavior in a multiband superconductors 60 . The 
corresponding fits are shown in Figs. 6 (a-c). Although 
the two-gap model contains two distinct gaps, the spe¬ 
cific heat value is calculated as the sum of contribu¬ 
tions, each one following the BCS-type temperature de¬ 
pendence, A(0) = 7iAi(0) + 72A 2 (0) 60 and the thermo¬ 
dynamic properties are obtained as the sum of the con¬ 
tributions from the individual bands, i.e., a\ = Ai /k B T c 
and a 2 = A 2 /k B T c . 

The estimated Ai(0)/fcBT c for the small gap for x 
= 0, 0.04, and 0.09 is 0.88T0.1, 1.9T0.2, and 2.2T0.2, 
while the large gap A 2 (0)/k B T c is found to be 2.2±0.2, 
2.5±0.2, and 2.35±0.2, for x = 0, 0.04, and 0.09 respec¬ 
tively. The calculated data and the relative weights are 
illustrated in red lines in Fig. 6 (a-c). The error bars repre¬ 
sents the width of the corresponding range of gap ampli¬ 
tudes obtained in the fit for both values of Ai(0)/fc_eT c 
and A 2 (0 )/k B T c . The results obtained in the present 
work are consistent with ones of the models considered 
in Ref. 15 . The ratio of the two gaps (A^Oj/A^O), is 
« 0.7 and 0.9 for x = 0.04, and 0.09, respectively) is 
comparable to the FeSeo. 43 Teo .57 case and it is notice¬ 
ably larger than in iron pnictide superconductors (be¬ 
tween 0.24 and 0.5) 37,38 . All of the fitting parameters 
are remarkably consistent with those obtained from the 
penetration depth measurements. They give a strong ev¬ 
idence for a two-gap SC at a Fe(Se,S) system. It has 
been theoretically demonstrated that in multiband su¬ 
perconductors if the ratio of two isotropic s-wave gaps 



FIG. 7: The field dependence of the mixed state quasipar¬ 
ticle contribution 7 (H) for H ||c for x=Q, 0.04, and 0.09. 
The dashed lines represent the phenomenological linear fits 
above H = IT. The upper inset presents the scaling of the 
data according to the s-wave scenario: C ca i -3 = [C(H) - 
C(0)]/T 3 vs. T/y/H. The lower inset shows the specific heat 
of FeSe0.96S0.04 plotted as C p /T vs. T 2 measured under var¬ 
ious magnetic fields up to 9 T in the low temperature region. 
The solid lines show a linear extrapolation of the data. 


A!(0)/A 2 (0) > 0.5, the field-induced low energy excita¬ 
tions would be less pronounced compared to a single¬ 
band s-wave symmetry 61 . With this respect, in the 
low field range j(H) would slowly increase with H (see 
Fig. 7). As mentioned above, our obtained ratio of the 
two gaps is higher compared to the critical value of 
0.5 suggested by the theory, which further confirms the 
multiband nature in Fe(Se,S). In addition, the obtained 
two gaps of both S-doped samples are consistent with the 
penetration depth results and larger than the BCS value 
in the weak-coupling regime. Overall, such a behavior 
confirms the strong coupling nodeless superconductivity 
in Fe(Se,S). 

Next we discuss the field dependence of specific heat, 
which is another independent, sensitive test of the gap 
structure. It has been well demonstrated that in the case 
of an isotropic s-wave superconductor, 7 (H) oc H be¬ 
cause the specific heat in the vortex state is dominated 
by the contribution from the localized QP in the vor¬ 
tex core. Recently, Storey et al. 62 pointed out that the 
number of Caroli-de Genned bound states increases lin¬ 
early with the field due to the linear increase in the num¬ 
ber of vortices entering the sample. On the other hand, 
for the line nodes 7 (H) oc H 0 - 5 , the QPs contributing 
to the density-of-states (DOS) come from regions away 
from the vortex core and close to the nodes and supercur¬ 
rents around the vortex core in the mixed state causing 
a Doppler shift of the QP excitation spectrum 63 . The 
temperature dependence of the low-T part of the specific 
heat data measured in various magnetic fields applied 
along the c axis is shown in the lower inset of Fig. 7. 
The data plotted as C p /T vs. T 2 fits to C p /T = + 
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/3T 2 , with 7 n and /? as electronic and lattice coefficients, 
respectively. It should be mentioned that the absence 
of the so-called at the linear-T term of the zero-field 
specific heat indicates high quality of the single crystals. 
Nearly perfect linear behavior without any magnetic im¬ 
purities have been observed in our samples (see Fig. 7). 
The applied magnetic field enhances the low-T specific 
heat, indicating the increase of the QP DOS at the Fermi 
level induced by a magnetic field. A linear extrapolation 
of the low-T data to zero temperature yields the field 
dependence of the field-induced contribution. The main 
panel of Fig. 7 presents the field dependence of the spe¬ 
cific heat coefficient. The dashed lines are linear fits for 
H || c above H = 1 T as anticipated for a case of nodeless 
SC gap. 

Further confirmation of the nodeless character of su¬ 
perconductivity in our investigated systems comes from 
the low temperature specific heat data of the finite- 
temperature region in the mixed state. In fact, the 
quasiparticle excitations in superconductors with differ¬ 
ent gap symmetries are obviously distinct. In s-wave 
superconductors, the inner-core states dominate quasi¬ 
particle excitations and a simple scaling law proposed by 
Liu et al. , holds in the case of possible s-wave gap in 
Sr 0 . 9 Lao.iCu 04 64 : 

C qp /T 3 « C core /T 3 = 7„/#c2(0) x (1 T/VH )" 2 (5) 

where Cqp and C cole are the specific heat of quasiparti¬ 
cles induced by the applied magnetic field and quasipar¬ 
ticles present inside from the Abrikosov vortex cores in 
the mixed state, respectively. The s-wave scenario of the 
scaling result of the field-induced term in the mixed state 
is presented in the upper inset of Fig. 7. All the data at 
different magnetic fields can be roughly scaled within the 
s-wave scenario in one line. Similar low temperature spe¬ 
cific heat studies have been already conducted on several 
Fe-based superconductors 65,66 . 

For the sake of comparison, in Table I we have sum¬ 
marized the superconducting parameters for the investi¬ 
gated samples extracted from this study. According to 
Fig. 1(d) and Table I, with increasing S content the criti¬ 
cal current density, which is a measure of the strength 
of the pinning force density and can be very conve¬ 
niently used to characterize the strength of disorder in 
the system, enhances, suggesting improved flux pinning 
in those samples. The absolute value of the penetration 
depth in the T —> 0 limit determined for FeSe 0 . 96 S 0 . 04 , 
FeSe 0 . 91 S 0 . 09 , and FeSe, yields A a6 (0) = 372(15), 433(15), 
and 446(15) nm, respectively. These values are some¬ 
what smaller than 560(20) nm found in Fe(Te,Se) 22 , but 
comparable to the FeSeo .85 and FeSe 18,19 . In Fig. 6(d- 
f), a kink structure is observed on the A~ 6 2 (T) curves. 
This kink in A A 2 (T) can be associated with the two- 
band supercondictivity as in the cases of Fe(Te,Se), Bao.6 
o. 4 Fe 2 As 2 , and MgB 2 22,67,68 . The upper critical field, 
H c 2 , for the S-doped FeSe sample increases with in¬ 
creased doping, which is mainly due to its enhanced T c . 


From 7 n values, we estimate the universal parameter 
AC e |/ 7 n T c , which is considerably higher than the pre¬ 
diction of the weak coupling BCS theory (A C e i/^ n T c = 
1.43). Taking into account the fact that the supercon¬ 
ducting transition is relatively sharp [see Fig. 6 (a-c)], a 
distribution in T c or the presence of impurity phases can¬ 
not explain the higher value of the universal parameter. 
We believe that the presence of strong coupling supercon¬ 
ductivity explain this higher values. Most remarkably, 
the specific-heat data allows for precise evaluation of SC 
volume fraction (V S c), i-e., V sc = (in ~ 1r)/ln, with 
7 r being the residual electronic specific-heat coefficient. 
Since our is almost absent (see lower inset in Fig. 7), 
VSC estimated from specific heat is in fair agreement 
with our magnetization data [see Fig. 1(a)]. The overall 
values of the investigated superconducting gap derived 
from specific heat is similar to the one obtained from 
the penetration depth. However, both large gap, A l, 
and smaller one, As, upon doping present a higher value 
than the weak-coupling BCS (1.76fcsT c ) gap value, which 
reflects a tendency for strong coupling effects. This is in¬ 
consistent with the theoretical constraints of the weakly 
coupled two-band superconductor model in which one 
gap must be larger than the BCS gap and one smaller 69 . 

Although, rather large single or multiple gap values 
were reported in Fe(Se,Te) from specific heat 23 , penetra¬ 
tion depth 22 , and ARPES 70 suggesting strong-coupling 
multiband superconductivity, the pairing symmetry in 
Fe(Se,Te) is still under debate. Additionally, two in¬ 
dependent reports of penetration depth measurements ' 1 
and scanning tunneling microscopy ' 2 in Fei+j / (Tei_ x Se a; ) 
have claimed the possibility of nodes in the SC gap. In¬ 
terestingly, near optimal doping FeSeo. 45 Teo. 55 , specific 
heat measurements demonstrate isotropic gap behavior 
under zero magnetic field but anisotropic/nodal gaps un¬ 
der magnetic field 23,39 . Our data show that Fe(Se,S) sys¬ 
tem belongs to the class of multiband superconductors, in 
the strong-coupling regime. Given the substantial diver¬ 
gency of the existing data on the gap values and the gap 
symmetry for FeSe-based superconductors, a combina¬ 
tion of several independent techniques rather than single 
technique is highly desirable. In the current paper we 
presented self-consistent data obtained from both lower 
critical field H c 1 , and specific heat measurements. We 
believe that other techniques such as fiSR , ARPES, and 
NMR are highly desirable to further confirm the multi¬ 
band structure in Fe(Se,S). 


IV. CONCLUSIONS 

In summary, using a AICI 3 /KC1 flux technique we have 
grown high quality single-crystals of FeSei-^S^ system 
(. x=0 , 0.04, 0.09, and 0.11) and studied their transport, 
magnetic and low temperature specific heat properties. 
We show that the introduction of S to FeSe enhances the 
upper critical field H c2 , critical current density J c , and 
the T c . The magnetic phase diagram has been studied 
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in the case of magnetic field applied along the c axis and 
ab plane for x = 0.04 and the resulting anisotropy was 
found to be around T = ~ 4. The temper¬ 

ature dependence of the penetration depth and C e i can 
be described neither within single band weak coupling 
BCS nor using the d -wave approach. Our results, (i) the 
T-dependencies of both penetration depth and specific 
heat, (ii) the kinky in A a b{T), (iii) the large specific heat 
SC gap values revealed from both probes, (iv) the linear 
field dependence of 7 , (vi) the large jump at T c , and (vii) 
the s-wave scaling of the low-T specific heat data in the 
mixed state, all indicate the presence of strong-coupling 
multiband and nodeless superconductivity in FeSei_ x Sa;. 
The field-induced change in the low-T specific heat shows 
a linear magnetic field dependence which is consistent 
with the s-wave symmetry of the order parameter. 
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